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Abstract 

The reaction pp — > pf(ir + 7r^TT°)p s has been studied at 450 GeV/c in an experiment designed 
to search for gluonic states. A spin analysis has been performed and the cIPt filter applied. 
, The analysis confirms the previous observation that all undisputed qq states are suppressed at 
small dPx- In addition, a clear difference is observed in the production mechanism for the rj 
H ' and u. 
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The WA102 collaboration has recently described the application of a kinematical filter on 
centrally produced events [0]. This filter, which has been proposed as a glueball-gg filter 0, 
was originally based on the idea of Double Pomeron Exchange (DPE) and hence it is interesting 
to observe the effect of this filter on channels that cannot be produced by DPE. This paper 
presents new results from the WA102 experiment on the centrally produced 7r + 7r~7r° final state 
in the reaction 

PP -»• Pf(n + ir'n°)p s (1) 

at 450 GeV/c. The subscripts f and s indicate the fastest and slowest particles in the laboratory 
respectively. Due to the G parity of the tt + ty^ty° system it cannot be produced by DPE. These 
data represent an increase by a factor of fifty in the previously published WA76 results on the 
centrally produced 7r + 7r~7T° final state 0. 

The data come from experiment WA102 which has been performed using the CERN Omega 
Spectrometer. The layout of the Omega Spectrometer used in this run is similar to that 
described in ref . [|J with the replacement of the OLG A calorimeter by GAMS 4000 || . 

Reaction (|1|) has been isolated from the sample of events having four outgoing charged 
tracks plus two 7s reconstructed in the electromagnetic calorimeter^ by first imposing the 
following cuts on the components of the missing momentum: |missing P x \ < 17.0 GeV/c, 
|missing P y \ < 0.16 GeV/c and |missing P z \ < 0.12 GeV/c, where the x axis is along the 
beam direction. A correlation between pulse-height and momentum obtained from a system of 
scintillation counters was used to ensure that the slow particle was a proton. 

The effective mass of the two 7s shows a clear 7r° signal (a = 17 MeV), which was selected by 
requiring 0.1 < m(77) < 0.17 GeV. The quantity A, defined as A = MM 2 (p f p s )~M 2 
where MM 2 (pfp s ) is the missing mass squared of the two outgoing protons, was then calculated 
for each event and a cut of |A| < 3.0 (GeV) 2 was used to select the 7r + 7r _ 7r° channel. Events 
containing a fast A ++ (1232) were removed if M(pfn + ) < 1.3 GeV, which left 1 323 136 centrally 
produced 7r + 7r~7r° events. After the selection of the 7r + 7T~7r° channel a kinematical fit was 
performed in order to apply overall energy and momentum balance. 

Fig. |l|a) shows the acceptance corrected n + n^n° effective mass spectrum renormalised to 
the total number of observed events. In addition to clear 7^(547) and cj(782) signals there is a 
broad enhancement at 1.2 GeV, some evidence for the a2(1320) and another broad enhancement 
at 1.67 GeV. 

A Dalitz plot analysis of the n + n^n° final state has been performed using Zemach tensors 
and a standard isobar model ||. The analysis has assumed p(770), a, / (980) and / 2 (1270) 
intermediate states and that only relative angular momenta up to 2 contribute. For the pn 
decays both 1 = and 1 = 1 final states have been considered. The a stands for the tity S-wave 
amplitude squared, and the parameterisation of Zou and Bugg has been used in this analysis. 

The geometrical acceptance of the apparatus has also been evaluated over the Dalitz plot 
of the 7r + 7r~7r° system in 40 MeV intervals between 0.8 and 2.0 GeV. In order to perform a 
spin parity analysis the log likelihood function, Cj = Z)j log Pj(i), is defined by combining the 

^^The showers associated with the impact of the charged tracks on the calorimeter have been removed from 
the event before the requirement of only two 7s was made. 
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probabilities of all events in 40 MeV n +r K 7r° mass bins from 0.84 to 2.0 GeV. In order to include 
more than one wave in the fit the incoherent sum of various event fractions aj is calculated: 

^ = Eiog(E%^W + ( 1 -E%)) ( 2 ) 

i \ 3 3 J 

where the term (1 — J2j a j) represents the phase space background which is a free parameter 
in each bin. The negative log likelihood function (— C) is then minimised using MINUIT [Q. 
Coherence between different J p states and between different isobar amplitudes of a given J p 
have been neglected in the fit. Different combinations of waves and isobars have been tried and 
insignificant contributions have been removed from the final fit. Using Monte Carlo simulations 
it has been found that the feed-through from one spin parity to another is negligible and that 
the peaks in the spin analysis cannot be produced by phase space or acceptance effects. The 
fit generates the phase space background as that part of the data not associated with a given 
wave on a bin by bin basis and one requirement is that this background is a smoothly varying 
function that shows no residual resonance structure. 

The result of the best fit is shown in fig. ^| for the total data sample. A fit using 

jPC = 

1++ p ± nT S-wave, J PC = 2++ pV D-wave, J PC = 2~+ p ± vr T P-wave, J PC = 2~+ / 2 (1270)7r° 
b-wave, J PC = 1 pn P-wave and phase space is found to be sufficient in order to describe 
the data. The phase space distribution resulting from the fit shows no remaining resonance 
structure. 

The J PC = 1 ++ p ± 7r =F S-wave shows a broad enhancement. It has been fitted using an 
expression of the form: 

dN _ m ai T ai (m) 



dm (m 2 — m? ai ) 2 + r„ (m)rn? ai 
where T ai (m) is a mass dependent width of the form 



^ (tti) 

F ai (m) = r ai m ai (4) 



where p 1+s {m) is the phase space for S-wave p ± (770)7r =F , including the effects of the mass 
dependent p ± (770) width and of interference between the two p ± (770) bands in the Dalitz 
plot |§ . The results of the fit are given in table [l] and the mass and width are consistent with 
the PDG values for the ai(1260) 0. 

The J PC = 2++ p^ix T D-wave shows a peak in the 02 (1320) region and has been fitted 
using a relativistic spin 2 Breit-Wigner. The results of the fit are given in table |l| The J PC = 
2~ + p^Ti T P-wave and J PC = 2 _+ / 2 (1270)7r° S-wave have similar distributions consistent 
with the 7T2(1670). The J PC = 2~ + p ± 7r =F P-wave has been fitted with a spin 1 relativistic 
Breit-Wigner and the J PC = 2 h /2(1270)7r° S-wave has been fitted with a spin relativistic 
Breit-Wigner. The masses and widths are similar and are given in table [l|. The J PC = 1 p ± 7r =F 
P-wave is dominated by a peak at the mass of the 0(1020) which has been fitted with a spin 1 
relativistic Breit-Wigner convoluted with a Gaussian to describe the experimental resolution. 

The 7r + 7r~7T° mass spectrum shown in fig. |I|a) has been fitted with the resonances described 
above together with a Gaussian (a = HMeV) to describe the 77, a Breit-Wigner convoluted with 
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a Gaussian (a = 18 MeV) to describe the uj and a background of the form a(m—m t h) b exp(— cm— 
dm 2 — em 3 ), where m is the 7r + 7r~7r° mass, m t h is the threshold mass and a,b,c,d,e are fit 
parameters. The fit is found to describe the data well and yields masses for the rj and uj given 
in table [[]. 

The branching ratio for the 7r 2 (1670) to pn and / 2 (1270)7r has been calculated taking into 
account the unseen decay modes of the / 2 (1270) and gives 

7r 2 (1670) -> fm nerr t n no+0.03 



tt 2 (1670) -> / 2 (1270) 



7T 



0.57 ±0.021^ (5) 



which is in good agreement with the PDG values of 0.55 ± 0.08 |L0]]. The systematic error 
quoted has three major sources. The first comes from the error in the branching ratio of 
/ 2 (1270) to 7T7T [HJ. The second is due to the determination of the number of events in each 



decay mode; this has been estimated by comparing the ratio calculated from the total data 
sample compared to the ratio calculated from the fits performed as a function of dPr (see 
below). Finally estimates of possible backgrounds below the signals have been made for both 
decay modes. 

In previous analyses of other channels it has been observed that when the centrally produced 
system has been analysed as a function of the parameter cIPt, which is the difference in the 
transverse momentum vectors of the two exchange particles [p]] , all the undisputed qq states are 
suppressed at small dPr- Therefore a study of the 7r + 7r~7r° system has been made as a function 
of dPx- The 7r + 7r~7r° mass spectrum is presented in fig. |l]b), c) and d) for dPx < 0.2 GeV, 
0.2 < dPr < 0.5 GeV and dPx > 0.5 GeV respectively. In order to determine the dependence 
of the 7] and the w as a function of dPx-, a fit has been made to these mass spectra using the 
parameters from the fit to the total mass spectrum. For the ai(1260), a 2 (1320), 7r 2 (1670) and 
0(1020) the spin analysis has been performed as a function of dPx and the results are shown 
in fig. [| The J PC = 2~ + p ± 7r =F P-wave and J PC = 2~ + / 2 (1270)7r° S-wave have been summed 
together to produce a single J PC = 2~ + wave. The waves have been fitted with the Breit- 
Wigners described above and the results are given in table |^. As can be seen all the resonances 
are suppressed at small dPx- However, it is interesting to note that they do not all have the 
same dPx behaviour. 

The azimuthal angle (<p) between the fast and slow proton is related to dPx by 

cost = 

where Pt is the transverse momentum of the central system and Pj^ 2 is the transverse momen- 
tum of the exchanged particle. From table |^ it can be seen that the 77 and the u> have a similar 
dPx dependence; however, they have a dramatically different cj) dependence. Fig. f|a) and b) 
show the angle (p for 0.5 < m(7r + 7r~7r° ) < 0.6 GeV (rj region) and for 0.735 < m(7r + 7r _ 7r° ) < 0.835 
GeV (u region) respectively. This suggests that the r\ and the cu have very different produc- 
tion mechanisms. Further evidence for this comes from studying the four momentum transfer 
squared from one of the proton vertices which is shown in fig. |]c) and d) for the rj and u) 
respectively. The \t\ distribution for the uj is more steep than that for the rj and extends to 
\t\ = unlike the r\ distribution which appears to turn over. 

In order to investigate the centre of mass energy dependence of rj and uj production, the ratio 
of the acceptance corrected number of rjs and u;s has been calculated using data from the WA76 



4 



experiment which was performed at an incident beam energy of 85 GeV/c (y/s = 12.7 GeV) || 
and from this current experiment (-^/i = 29 GeV). 

^\ = 0.20 ± 0.02 (v^ = 12.7 GeV) = 0.09 ± 0.01 (x/s = 29 GeV) (6) 
a(u) 

which again indicates a different production mechanism for the two resonances. 

In summary, a study of the centrally produced n + n~TT system shows that the most promi- 
nent resonance signals are due to the r/, to, 0(1020), a! (1260), a 2 (1320) and 7r 2 (1670). All these 
states are found to be suppressed at small dP?. In addition, a clear difference is observed in 
the production characteristics of the r] compared to the u. 
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Tables 



Table 1: Parameters of resonances in the fit to the tt + ty n mass spectrum and waves. 





Mass (MeV) 


Width (MeV) 


Observed 
decav mode 


I{J PC ) 


Number of 
events 




549 ± 1 






0(0 _ +) 


21590±644 


cj(782) 


785 ± 2 


8.4 (fixed) 




0(1") 


248511±2118 


0(1020) 


1021 ± 4 


4 (fixed) 


pTY 


0(1") 


9053±490 


Oi(1260) 


1240 ± 10 


400 ± 35 


pTY 


1(1++) 


327957±2678 


a 2 (1320) 


1317 ± 3 


120 ± 10 


pTY 


1(2++) 


35616±648 


vr 2 (1670) 


1669 ± 4 

1670 ± 4 


268 ± 15 
256 ± 15 


pn 
/ 2 (1270)tt 


1(2-+) 
1(2-+) 


23535±598 
23692±549 
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Table 2: Resonance production as a function of dPx expressed as a percentage of its total 
contribution. 





dP T <0.2 GeV 


0.2<dP T <0.5 GeV 


dP T >0.5 GeV 


V 


5.3 ± 0.5 ± 0.5 


42 ± 1 ±1 


52 ± 1 ± 1 




13 ± 1 ± 1 


44 ± 1 ± 1 


43 ± 1 ± 1 


0(1020) 


8 ± 2 ± 2 


47 ± 3 ± 3 


45 ± 2 ± 4 


ai(1260) 


17 ± 1 ± 3 


53 ± 3 ± 2 


29 ± 2 ± 2 


a 2 (1320) 


4 ± 2 ± 3 


37 ± 2 ± 1 


59 ± 3 ± 3 


tt 2 (1670) 


15 ± 1 ± 1 


51 ± 2 ± 3 


33 ± 1 ± 3 
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Figures 



Figure 1: The ir + n 7r° effective mass spectrum a) for the total data with fit using 3 Breit- 
Wigners b) for dP T < 0.2 GeV, c) for 0.2 < dP T < 0.5 GeV and d) for dP T > 0.5 GeV. 



Figure 2: Results of the spin parity analysis. The superimposed curves are the resonance 
contributions coming from the fits described in the text. 



Figure 3: Results of the spin parity analysis as a function of dPr- 



Figure 4: The azimuthal angle (</>) between the two outgoing protons for the a) rj and b) uj. 
The four momentum transfer squared (|t|) from one of the proton vertices for the c) rj and d) 

UJ. 
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